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ABSTRACT 

We present a study of correlations between the elemental abundances and galaxy morphologies of 
91 blue compact galaxies (BCGs) at z=0. 20-0.35 with Sloan Digital Sky Survey (SDSS) DR7 data. 
We classify the morphologies of the galaxies as either 'disturbed' or 'undisturbed', by visual inspec- 
tion of the SDSS images, and using the Gini coefficient and M2o- We derive oxygen and nitrogen 
abundances using the T e method. We find that a substantial fraction of BCGs with disturbed mor- 
phologies, indicative of merger remnants, show relatively high N/O and low O/H abundance ratios. 
The majority of the disturbed BCGs exhibit higher N/O values at a given O/H value compared to 
the morphologically undisturbed galaxies, implying more efficient nitrogen enrichment in disturbed 
BCGs. We detect Wolf-Rayet (WR) features in only a handful of the disturbed BCGs, which appears 
to contradict the idea that WR stars are responsible for high nitrogen abundance. Combining these 
results with Galaxy Evolution Explorer (GALEX) GR6 ultraviolet (UV) data, we find that the ma- 
jority of the disturbed BCGs show systematically lower values of the Ha to near-UV star formation 
rate ratio. The equivalent width of the H/3 emission line is also systematically lower in the disturbed 
BCGs. Based on these results, we infer that disturbed BCGs have undergone star formation over 
relatively longer time scales, resulting in a more continuous enrichment of nitrogen. We suggest that 
this correlation between morphology and chemical abundances in BCGs is due to a difference in their 
recent star formation histories. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: starburst — galaxies: star 
formation 



1. INTRODUCTION 

Blue compact galaxies (BCGs) are compact, star- 
bursting galaxies which are thought to be the least chemi- 
cally evolved ones in the local Universe, and are therefore 
considered as analogs of high-redshift star-forming galax- 
ies. They allow important constraints to be placed on 
massive star formation (SF) and the evolution of galaxies 
in very low metallicity environments (Thuan & Martin 
1981; Izotov & Thuan. 1999). While a number of studies 
have probed the formation and evolution of BCGs, merg- 
ers and interactions between various types of galaxy have 
been suggested as mechanisms to stimulate enhancement 
of SF activity in BCGs (e.g., Larson & Tinsley 1978; Tay- 
lor et al. 1995, Mendez & Esteban 2000; Pustilnik et al. 
2001a, b; Ostlin et al. 2001; Sung et al. 2002; Johnson 
et al. 2004; Ekta & Chengalur 2010a). Recently, based 
on numerical simulations, Bekki (2008) showed that blue 
compact dwarf (BCD) galaxies can be formed by mergers 
and interactions between gas-rich dwarfs. In his simula- 
tions, mergers are able to trigger a central starburst to 
form young, massive cores, resembling the morphologies 
of currently known BCD galaxies. 

The gas phase elemental abundances of BCGs is a 
key physical parameter which traces their star forma- 
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tion history (SFH) and chemical evolution. The N/O 
abundance ratio as a function of the O/H ratio is con- 
sidered to be an important observational constraint for 
probing the chemical evolution of galaxies (Edmunds & 
Pagel 1978; Torres-Peimbert, Peimbert, & Fierro 1989; 
Vila-Costas & Edmunds 1993; Izotov et al. 2006). There 
has been observational evidence that the oxygen and ni- 
trogen abundances may be profoundly linked to the local 
environments of galaxies. Many studies have suggested 
that mergers and tidal interactions can induce a low O /H 
ratio in the galaxies' central region owing to inflows of 
metal-poor gas from the outer part of the galaxy (Rupke 
et al. 2008, 2010; Peeples et al. 2009; Ekta & Chen- 
galur 2010b; Kewley et al. 2010; Montuori et al. 2010). 
While the N/O ratio was found to be fairly constant for 
low-metallicity BCGs (Izotov & Thuan 1999), a nitro- 
gen overabundance was detected in some BCGs. Fur- 
thermore, some studies suggested that the unexpectedly 
high N/O ratios of BCG could be associated with the 
merger process and hence with the related powerful star- 
burst (e.g., Pustilnik et al. 2004). Recently, based on a 
sample of star-forming luminous compact galaxies at z < 
0.63, Izotov, Guseva, & Thuan (2011) found a significant 
spread in N/O values, with a hint of enhanced nitrogen 
abundance. Furthermore, in their sample, a considerable 
fraction of the galaxies show signs of having disturbed 
morphologies, suggesting recent interactions. 

All of these previous results suggest that the effects of 
mergers and interactions provide a remarkable opportu- 
nity for clarifying the origin of the scatter in the N/O vs. 
O/H diagram. However, despite the importance of ele- 
mental abundances, detailed studies of their relationship 
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with the morphology and the SFH of BCGs have been 
limited by small samples. To improve this situation, in 
this letter, we examine a large and homogeneous sample 
of BCGs at intermediate redshifts around z~0.3. Our 
goal is to assess whether the BCGs are discriminated by 
features of their elemental abundances, linking their SF 
processes to the merger and interaction of galaxies. 

2. SAMPLE SELECTION AND ANALYSIS 

We select a sample of 91 highly excited BCGs in the 
rcdshift range of 0.20 < z < 0.35 using criteria that are 
based on both spectra and images from the Sloan Dig- 
ital Sky Survey (SDSS) Data Release 7 (DR7) (Abaza- 
jian et al. 2009). Only galaxies with a well-detected 
[OIIIJA4363 emission line, with S/N higher than 5, are 
selected. This criterion allows a direct determination of 
the electron temperature and abundances. We select 
only galaxies which occupy the region of star-forming 
galaxies with the highest excitation and which are well 
separated from active galactic nuclei in the BPT diagram 
(Baldwin, Phillips, & Terlevich 1981). In addition, only 
BCG galaxies with relatively compact morphologies and 
blue colors in their SDSS images are considered. Our 
BCGs span a luminosity range of -22 < M r < -19 and 
color ranges of 0.1 < g — r < 1.2 and -2.0 < r — i < 
0.2. According to the photometric selection criteria used 
by Cardamone et al. (2009), about 59 % (54 of 91) of 
our BCGs are classified as green pea galaxies. Our sam- 
ple also includes 14 BCGs which show Wolf-Rayet (WR) 
features based on the detection of a blue bump around 
4650A and/or the He II emission line at 4686A (Conti 
1991; Pustilnik et al. 2004; Brinchmann, Kunth, & Dur- 
ret 2008). 

By visual inspection of the SDSS g, r, i combined color 
images, we conducted a systematic morphological analy- 
sis in which we classified each galaxy as either disturbed 
or undisturbed. We classified 19 galaxies as disturbed, 
on account of their asymmetric and distorted outer light 
distribution and/or cometary appearance, indicative of 
merger remnants or recent disturbance by tidal inter- 
action. The 72 BCGs classified as undisturbed show 
compact morphologies of relatively round and smooth 
shape. Example SDSS images for undisturbed and dis- 
turbed BCGs are shown in Figure 1. For the four undis- 
turbed BCGs in our sample which also have an optical 
HST or SDSS Stripe 82 image, we have tested the impact 
of image quality on our morphological classification (see 
inset images of Fig. la, b for two objects). None of these 
higher-quality, ancillary images show clear evidence of 
mergers or interactions, although one does show an elon- 
gated shape with a hint of two separate SF regions (see 
Fig. la). Note that the disturbed BCGs are system- 
atically biased towards higher luminosities (-22 < M r < 
-21), while undisturbed BCGs span the whole luminosity 
range of our sample. 

In addition to our visual classification, we have de- 
rived the Gini coefficient and M20 from the SDSS r- 
band images, using the method of Lotz, Primack, & 
Madau (2004). In G-M20 parameter space (Figure 2), 
the disturbed and undisturbed BCGs show clearly dif- 
ferent distributions. In particular, the disturbed BCGs 
show higher G and lower M20. In other words, the dis- 
turbed BCGs show steeper, and thus more centrally- 
concentrated intensity distributions than the undis- 



turbed BCGs, which is the result of distorted low surface- 
brightness regions in the outskirts of the disturbed galax- 
ies. This result supports our visual morphological classi- 
fication. 

From the Galaxy Evolution Explorer (GALEX) GR6 
data, 35 BCGs are detected in the far-ultraviolet 
(FUV; 1350-17501) and the near-ultraviolet (NUV; 
1750-2800A) bands. The FUV and NUV fluxes are taken 
from the GALEX pipeline data. Based on the redden- 
ing law of Cardelli et al. (1989), galactic extinction was 
corrected using the E(B-V) value from the NED. The in- 
ternal reddening coefficient C(H/3) was derived from the 
SDSS spectrum and its absorption was determined from 
A(FUV) = 8.19 x E(B-V) = 5.72/C(H/3) and A(NUV) 
= 9.21 x E(B-V) = 6.44/C(H/3). 

The global fluxes of the emission lines from each galaxy 
were systematically measured by fitting a Gaussian pro- 
file, using the IRAF splot task. We applied a foreground 
reddening correction to the line intensities following Lee 
& Skillman (2004). Using the traditional Balmer decre- 
ment method (Izotov et al. 1994), the observed line in- 
tensities were also corrected for both internal redden- 
ing and underlying stellar absorption. The gas-phase 
elemental abundances of O and N were derived by the 
direct method, using the STSDAS nebular /ionic task. 
Using the ionization correction factors from Izotov et al. 
(2006), we corrected for unobserved ionization stages and 
then estimated the total elemental abundances. 

We derive the Ha, FUV, and NUV star formation rates 
(SFRs) of our BCGs using the relations given by Ken- 
nicutt et al. (1998). We have corrected the Ha flux to 
take into account light-loss from the 3 arcsec diameter 
of the SDSS spectroscopic aperture. The total Ha flux 
was estimated using the ratio between Petrosian and 3 
arcsec aperture fluxes in the i-band surface photometry. 
The SDSS aperture collects on average 80% of the to- 
tal Ha flux of our galaxies, owing to the compactness of 
the galaxies (i-band Petrosian radii are in the range 0.9" 
-3.3"). In Table 1, we summarize the median values of 
basic parameters derived from our analysis, for disturbed 
and undisturbed BCGs. 

In addition, a sample of 104 local BCGs, with a sim- 
ilar luminosity range to our BCG sample, was compiled 
from the literature (Izotov & Thuan 1999; Izotov et al. 
2006) in order to compare the properties of BCGs at low 
and intermediate redshifts. Using SDSS DR7 images, we 
classified the morphologies of the local BCGs into dis- 
turbed and undisturbed galaxies as described above. The 
gas-phase abundances of local BCGs are also adopted 
from Izotov & Thuan (1999) and Izotov et al. (2006). 
For comparison with our disturbed BCGs, we also se- 
lected five well-known examples of local BCGs show- 
ing clear merger features (HS 0837+4717, III Zw 107, 
UM420, SBS 0926+606, and IRAS 08208+2816; Pustil- 
nilk et al. 2004; Lopez-Sanchez 2010). All of these local 
merger BCGs are detected by GALEX, from which we 
derive their FUV and NUV SFR. The Ha SFR of all but 
one (HS 0837+4717) are available from Kennicutt et al. 
(2008). 

3. RESULTS 
3.1. N/O and O/H ratios 
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In Figure 3, we show N/O vs. O/H for the disturbed 
(red filled circles) and undisturbed (black filled circles) 
BCGs. For comparison, we also show the distributions 
of disturbed (solid contours) and undisturbed (dashed 
contours) local BCGs. The BCGs show metallicities in 
the range of 7.6 < 12+log O/H < 8.5 with a median 
value of 8.13 ± 0.37. This corresponds to <~ 27% of the 
solar value assuming 12+log(0/H)© = 8.69 (Allende Pri- 
eto et al. 2001), indicating that these systems are gen- 
uine metal-poor galaxies. In the case of the galaxies with 
12+log(0/H) > 8.0, the overall distribution of our BCGs 
shows no prominent difference to that of local BCGs, in- 
dicating no segregation of elemental abundances of BCGs 
at z < 0.35 (see also Izotov et al. 2011). On the other 
hand, it is worth noting that our BCGs with low O/H 
values (12 + log O/H < 8.0) appear to have system- 
atically higher N/O values compared to those showing 
relatively higher O/H values. 

Regarding the distribution of elemental abundances 
as a function of morphology, one remarkable result is 
that the disturbed BCGs exhibit a clearly different dis- 
tribution to that of the undisturbed ones. Most of the 
disturbed BCGs display systematically higher N/O val- 
ues at a given O/H value compared to the undisturbed 
ones. The mean difference of N/O values between two 
subsamples is very pronounced (about 0.2 dex, see Ta- 
ble 1 and histograms in the right panel of Fig. 3). A 
K-S test on the N/O distributions between disturbed 
and undisturbed BCGs yields a probability of 0.04 % 
for the same parent distribution of the two subsam- 
ples. At the same time, a large fraction of the disturbed 
BCGs is also systematically biased towards low O/H val- 
ues (12+log(0/H) < 8.0). In line with previous results 
(Rupke et al. 2008, 2010; Peeples et al. 2009; Ekta 
& Chengalur 2010b; Kewley et al. 2010; Montuori et al. 
2010), the low O/H values of some of the disturbed BCGs 
support the idea that these galaxies are associated with 
recent mergers or interactions with neighboring galax- 
ies. When we consider only the disturbed BCGs, there 
appears to be an anticorrelation between the N/O and 
O/H ratios such that disturbed BCGs with low O/H val- 
ues show systematically high N/O values. 

Note that the five local disturbed BCGs (blue squares 
in Fig. 3) exhibit a similar distribution to our disturbed 
BCGs, and most of them also have high N/O values. 
The most extreme object is HS 0837+4717, which is a 
very metal-deficient BCG (12+log(0/H)=7.64) showing 
a significant N excess (log(N/0)= -0.83: Pustilnilk et al. 
2004). 

3.2. Wolf-Rayet Features 

The bias towards high N/O values in morphologi- 
cally disturbed BCGs necessitates mechanisms for the 
enhanced production of nitrogen in these galaxies. One 
possibility is that WR stars are formed in starbursts trig- 
gered by mergers. Since merging of two galaxies provides 
a subsequent intense starburst, this could result in a sig- 
nificant increase of the relative fraction of massive young 
stars. Four of the five local merger BCGs have previously 
been shown to have clear WR features (blue squares with 
yellow cross in Fig. 3; Pustilnilk et al. 2004; Lopez- 
Sanchez 2010, but see James, Tsamis, & Barlow 2010 for 
non-detection of WR features in UM420 from their IFU 
observations). However, WR features are detectable for 



only a handful of our disturbed BCGs (four red filled 
circles with yellow cross in Fig. 3). Furthermore, WR 
features are not preferentially found in the sample of dis- 
turbed BCGs with high N/O values, regardless of mor- 
phology. 

On the other hand, three of the four disturbed BCGs 
with WR features appear to show morphologies more 
typical of an on-going merger or approaching post-merger 
(i.e., contact between two separate galaxies or distinct 
tidal tail, see Fig. Id, c, f), although it is somewhat dif- 
ficult to properly assess detailed morphologies with just 
the SDSS images. In the case of the local BCGs with 
WR features, they also show similar, on-going merger 
morphologies (e.g., two separate galaxies for III Zw 107 
and UM420 as well as a distinct, long tidal tail for IRAS 
08208+2816). This suggests that the appearance of WR 
features in disturbed BCGs might be related to the de- 
gree of disturbance in their morphologies, which itself is 
an indicator of merger-status and, therefore, starburst 
activity. 

The detection and discrimination of WR features in 
the spectrum of a star burst galaxy can be hampered by 
several effects, including the quality of the spectrum, or 
the location and size of the spectroscopic aperture (see 
Lopez-Sanchez & Esteban 2008 and references therein). 
In the case of large apertures, weak WR features can be 
strongly diluted by the continuum flux (Lopez-Sanchez 
& Esteban 2008; James et al. 2010). For this reason, 
we cannot rule out the existence of WR populations in 
the disturbed BCGs, even though they are not clearly 
detected in the SDSS spectra. 

3.3. Star Formation Histories 

The starburst duration directly affects the chemical 
yield and the stellar age of a galaxy. Long bursts of 
SF are able to produce high N/O values, without nec- 
essarily showing distinct WR features. While WR stars 
are widely separated in time, resulting in a weak or no 
WR features in the spectrum, the interstellar medium 
(ISM) in BCGs is cumulatively enriched over time by 
massive stars. Since the disturbed BCGs in our sample 
are systematically biased towards higher masses (-22 < 
M r < -21) (see also Telles, Melnick, & Terlevich, 1997), 
they can retain a large fraction of newly produced met- 
als. The relatively long duration of the SF for BCGs 
that are formed by the merger of very gas-rich dwarfs is 
supported by numerical simulations (200 - 300 Myr, see 
Fig. 4 of Bc'kki 2008). 

Different indicators of the SFR are sensitive to stars of 
different masses and timescales. The Ha emission traces 
the present SF activity from massive, early type stars 
with timescales of a few tens of Myr. On the other hand, 
the UV flux traces less massive stars with timescales of 
a few hundred Myr (Kennicutt 1998; Iglesias-Paramo et 
al. 2004) . The Ha luminosity declines more rapidly than 
the UV continuum luminosity, which leads to a decrease 
in the Ha to UV SFR ratio with increasing starburst 
age (Sanchez-Gil et al. 2011). Thus, the Ha to UV 
SFR ratio carries important information about the SFH 
and the duration of the SF episode of a galaxy. In Fig- 
ure 4 (left panels), there is a striking difference between 
the disturbed (red filled circles) and undisturbed BCGs 
(black filled circles) , in terms of their Ha to NUV SFR 
ratios; the disturbed BCGs show SFR ratios that are 
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about a factor of 5 lower than the undisturbed BCGs 
(see also Table 1). Interestingly, the nearby BCGs that 
show merger-related features (blue squares) are also con- 
sistent with this result. This supports the idea that dis- 
turbed BCGs have undergone SF over relatively longer 
time scales than the undisturbed BCGs. 

The H/3 emission line equivalent width, EW(H/3), has 
also been suggested as an indicator of the age of the lat- 
est burst of SF in a galaxy (Stasinska & Izotov 2003; 
Izotov et al. 2011). It is clearly seen that the majority 
of the disturbed BCGs also show systematically lower 
EW(H/3) than the undisturbed ones (about a factor of 
3), implying a significant underlying, non- ionizing older 
stellar population which contributes to the stellar contin- 
uum (see right panels of Fig. 4 and Table 1). This result 
provides further support for the idea that disturbed and 
undisturbed BCGs have different SFH. 



4. DISCUSSION AND CONCLUSIONS 

We have established that the chemical abundances 
within BCGs at 0.20 < z < 0.35 are closely correlated 
with galaxy morphology. The BCGs in our sample are 
segregated by their morphology in the N/O vs. O/H di- 
agram. While there have been several observational and 
theoretical suggestions that galaxy mergers or interac- 
tions result in a decrease in O/H, our study has shown 
that a significant fraction of morphologically disturbed 
BCGs also exhibits nitrogen enhancement. Disturbed 
BCGs with an elevated nitrogen abundance might have 
experienced prolonged SFHs through dynamical events 
(e.g., merger or interaction) between two galaxies, and 
consequently the ISM is continuously enriched by mas- 
sive stars. Considering the high fraction of green pea 
galaxies in our BCGs, our results are consistent with re- 
cent studies of abundance patterns in green pea galax- 
ies, which have revealed systematically higher N/O ra- 
tios compared to most star forming galaxies at a given 
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O abundance, possibly resulting from recent interaction- 
induced inflows of gas (Amorin et al. 2010, 2012). Inflow 
of very low metallicity gas may dilute O/H, without sig- 
nificantly affecting the N/O ratio. This might explain 
the simultaneously high N/O and relatively low O/H ra- 
tios of some of the BCGs, due to the horizontal offset 
towards lower O/H values in the N/O vs. O/H diagram 
(Koppen & Hensler 2005). On the other hand, it is im- 
portant to point out that our results are not affected by 
the difference in luminosity of the disturbed and undis- 
turbed subsamples. If we instead consider only the BCGs 
in the luminosity range -22 < M r < -21, the results of 
our analysis are not significantly changed. 

Our results pertain to the presence of the strongest 
optical morphological signatures of merging or interac- 
tion. However, due to the limited spatial resolution (i.e., 
8 kpc resolution for the 1.4 arcsec SDSS PSF at z - 0.3) 
and depth of the SDSS images used, we cannot rule out 
that we have missed faint merging/interaction morpho- 
logical features in some of the apparently undisturbed 
BCGs. Based on the trends identified in our sample, ob- 
jects with low level, small scale merging or interaction 
features may appear as morphologically undisturbed in 
the SDSS images but with a high N/O, low Ha to UV 
SFR ratio, and low EW(H/3). Further deep, high resolu- 
tion imaging observations will be necessary for clarifying 
this issue. 
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J0820+5050 z=0.2l J22I9+000I z=0.23 




Fig. 1. — Example of 12.6" X 12.6" SDSS g, r, i composite color images of morphologically undisturbed and disturbed BCGs: (a, b) 
undisturbed BCGs with round and compact shapes. Inset images are those from HST/WFPC2 or SDSS Stripe 82 field; (c) disturbed BCG 
without WR features in its spectrum; (d, e, f) disturbed BCGs with WR features in their spectra. In each panel, the SDSS name and 
redshift of the galaxy is indicated. 
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Fig. 2. — G vs. M20 diagram for disturbed (red circles and contours) and undisturbed (black circles and contours) BCGs. Error bars 
indicate the mean values and standard deviations of G and M20 for these two subsamples, while the histograms show their distributions. 
The mean values of G and M20 are [0.60, -1.79] and [0.56, -1.58] for disturbed and undisturbed BCGs, respectively. 
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Fig. 3. — N/Ovs. O/H diagram of BCGs. Red and black filled circles are disturbed and undisturbed BCGs, respectively. The distributions 
of the disturbed and undisturbed local BCGs from the literature (Izotov &; Thuan 1999; Izotov et al. 2006) are denoted by solid and dashed 
contours, respectively. Blue squares represent the local merging BCGs. BCGs with WR features are represented by crosses. The histograms 
show the distributions of N/O and O/H for the disturbed and the undisturbed BCGs. 
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Fig. 4. — (a, c) Elemental abundance vs. the Ho to NUV SFR ratio of the BCGs. (b, d) Elemental abundance vs. H/3 emission line 
equivalent width, EW(H/3), of the BCGs. The dashed line in each panel presents the median value of the undisturbed BCGs. Symbols are 
the same as in Figure 3. The histograms show the distributions of He? to NUV SFR ratio and EW(H/3) for the disturbed and undisturbed 
BCGs. 
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TABLE 1 

Average characteristics of the disturbed and undisturbed BGCs 



Sample 


12+log(0/H) 


log(N/0) 


SFR Ha 


SFR]v(7V 


log [SFRWSFR W y] 


EW(H/3) 








[M /yr] 


[M /yr] 




[A] 


Disturbed BCG 


8.11 


-1.15 


10.91 


50.43 


-0.64 


48 


Undisturbed BCG 


8.15 


-1.35 


6.99 


8.82 


0.11 


136 



